Superconducting nanofilms: molecule-like pairing induced by quantum confinement 
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Quantum confinement of the perpendicular motion of electrons in single-crystalline metallic su- 
perconducting nanofilms splits the conduction band into a series of single-electron subbands. A 
distinctive feature of such a nanoscale multi-band superconductor is that the energetic position of 
each subband can vary significantly with changing nanofilm thickness, substrate material, protection 
cover and other details of the fabrication process. It can occur that the bottom of one of the avail- 
able subbands is situated in the vicinity of the Fermi level. We demonstrate that the character of 
the superconducting pairing in such a subband changes dramatically and exhibits a clear molecule- 
like trend, which is very similar to the well-known crossover from the Bardeen-Cooper-Schrieffer 
regime to Bose-Einstein condensation (BCS-BEC) observed in trapped ultracold fermions. For Pb 
nanofilms with thickness of 4 and 5 monolayers (ML) this will lead to a spectacular scenario: up to 
half of all the Cooper pairs nearly collapse, shrinking in the lateral size (parallel to the nanofilm) 
down to a few nanometers. As a result, the superconducting condensate will be a coherent mixture 
of almost molecule-like fermionic pairs with ordinary, extended Cooper pairs. 

PACS numbers: 74.78.Na 



I. INTRODUCTION 

Advances in self-assembly and nanofabrication resulted 
recently in the thinnest superconductors, i.e., single- 
crystalline metallic nanofilms with atomically uniform 
thickness down to a few monolayers, see, e.g., Refs. Q]-[5|. 
They are characterized by an extraordinary purity and, 
so, the only issue of disorder in the case of interest is due 
to the interface between the substrate and the nanofilm. 
Such an interface has an effect on the reflection of elec- 
tron waves (it is not perfectly specular). Nevertheless, 
there are clear signatures of the formation of discrete lev- 
els for the perpendicular electron motion [i.e., quantum- 
well states (QWS)] in tunneling spectral Pairing cor- 
relations appear to be very robust in single-crystalline 
nanofilms so that even Pb and In films with thickness of 
1 ML were found not to exhibit any considerable signs 
of degradation of superconductivity due to fluctuations. 4 
Experimental results for the temperature dependence of 
the pairing gap show almost no deviations from the or- 
dinary BCS picture, see, e.g., Refs. tJJMJ- However, it is 
worth noting that, for some samples, there appear sig- 
natures of suppression of density of states (DOS) in the 
tunneling spectra at low voltage above T c , which resem- 
bles, to some extent, the pseudogap physics of high-T c 
superconductors 4^ 

Another recent experimental breakthrough concerns 
the iron-pnictides (for a review see, e.g., Ref. 0), i.e., 
a new generation of multi-band superconductors spec- 
ified by the interplay of different band-dependent su- 
perconducting condensates. This renewed the interest 
in multi-band superconductivity. In particular, one of 
the most intriguing points is a competition of the char- 
acteristic lengths of the different condensates^— and a 
possible contribution of such a competition to new phe- 
nomena, including recently observed and strongly de- 



bated unconventional patterns of the vortex distribution 
in magnesium-diboride. 9 

The question arises if there is any relation between 
bulk multi-band superconductors and high-quality super- 
conducting nanofilms? In fact, they are similar in many 
important respects. As already mentioned above, tunnel- 
ing experiments demonstrate that the conduction band 
in metallic nanofilms splits up into a series of subbands 
in such a way that the bottom of each subband is at 
the energy position of the corresponding QWSji^ For 
example, Pb nanofilms with thickness 4-5 ML are two- 
band(subband) superconductors i^ In those metallic su- 
perconducting nanofilms the energetic position of each 
subband with respect to the Fermi level Ep changes with 
thickness, i.e., the discrete perpendicular levels scale as 
1/cP, with d the thickness of the nanofilm while the ac- 
companying variation of Ep is almost insignificant. 22 Ad- 
ditional reasons for shifting in energy of the subband bot- 
toms as measured from Ep are due to fabrication circum- 
stances: the presence or absence of a protective cover, 
an effect of the wetting layer, use of different substrates 
etc., which can significantly change the behavior of the 
single-electron wave functions at the interface (and, so, 
the QWS energy). Thus, the relevant microscopic param- 
eters of the different subbands, e.g., DOS and the Fermi 
velocity (in the parallel direction), are not fixed and can 
significantly vary from sample to sample even for the 
same thickness. It is possible to expect that these param- 
eters can, in principle, be tailored in future experiments. 
Therefore, it is timely to investigate in more detail the 
properties of the superconducting state in metallic single- 
crystalline nanofilms in the case when the bottom of one 
of the relevant single-electron subbands approaches the 
Fermi level. Below we show that the pairing in such 
a subband exhibits a clear molecule-like trend, which is 
very similar to the BCS-BEC crossover investigated at 



length in cold atomic gases, see, e.g., Ref. Il2|. In par- 
ticular, we predict that striking results can be obtained 
for Pb nanofilms with thicknesses 4-5 ML, where up to 
half of all the Cooper pairs nearly collapse, i.e., shrink- 
ing in the lateral size (parallel to the nanofUm) down to 
a few nanometers. Our finding significantly compliments 
the recent first observation of the BCS-BEC crossover in 
a solid-state material, i.e., in one of the available sub- 
bands in a multiband superconducting iron-chalcogenide 
FeSe x Tei_ x , see Ref. [H 

The present paper is organized as follows. In Sec. UU 
we address the criterion of the BCS-BEC crossover-like 
behavior in a given single-electron subband. We show 
that a condensed pair of electrons significantly shrinks in 
its lateral size (parallel to the nanofilm) when the ratio 
of the parallel kinetic energy of electrons to the abso- 
lute value of the potential energy becomes smaller than 
one, i.e., when the subband bottom approaches the Fermi 
level. In Sec. IIIII we investigate 4 and 5 ML thick lead 
nanofilms, where only two single-electron subbands are 
occupied, i.e., we have a coherent mixture of two con- 
densates, and one of them is specified by an extremely 
small characteristic length. Our conclusions are given in 

sec. inn 



II. BCS-TO-BEC CROSSOVER DRIVEN BY 
QUANTUM-SIZE EFFECTS 

It is well-known that the governing parameter for the 
BCS-BEC crossover in the system with fermionic pair- 
ing correlations is the ratio of the relevant kinetic en- 
ergy K to the absolute value of the potential energy U . 
When K/\U\ > 1, we are in the BCS regime with loosely 
bound and rather extensive Cooper pairs (the BCS limit 
corresponds to K/\U\ — > oo). On the opposite side of 
the crossover, when K/\U\ <C 1, molecule-like bound 
pairs appear with the BEC limit K/\U\ — > 0. Here it 
is worth noting that an extensive analysis of the behav- 
ior of the kinetic and potential energy through the BCS- 
BEC crossover for the attractive Hubbard modet^ gave 
K/\U\ « 2.2 at the BCS side, while approaching the 
BEC regime it resulted in K/\U\ « 0.03. Interestingly, 
the right balance between the strength of the kinetic and 
interaction energy has been considered as an essential 
feature of high-T c superconductivity. 15 

Keeping in mind the criterion based on the ratio 
K/\U\, let us examine what happens with a Cooper pair 
in a given subband when changing the position of its bot- 
tom with respect to Ep in single-crystalline nanofilms. 
If the bottom of the subband is situated far below the 
Fermi level, the mean kinetic energy of electrons in this 
subband comes mostly from the parallel motion, i.e., 
K ± < #|| - Ep. In this case K\\/\U\ ~ E F /\U\ > 
1. However, when the subband bottom (or the corre- 
sponding QWS) approaches the Fermi surface, we have 
K\_ — > Ep and, so, K\\ becomes much smaller than Ep. 
The parallel motion of electrons is reduced, and such 



a redistribution of the kinetic energy between the per- 
pendicular and parallel spatial degrees of freedom leads 
to a significant decrease in the ratio i^||/|t/|. Further- 
more, when the subband bottom goes above the Fermi 
level (such a subband still makes a contribution unless 
its bottom is above Ep + huo, with ujp, the Dcbye fre- 
quency), K\\ — > and, hence, J^n/IC/j becomes extremely 
small, i.e., A"||/|[/| <c 1. Therefore, when assuming that 
the ratio K\\/\U\ controls the lateral size of a condensed 
fermionic pair associated with the corresponding sub- 
band in the nanofilm, we expect from the above argu- 
ments a significant reduction in the lateral size of the 
Cooper pairs. 

To go into more detail, let us consider a superconduct- 
ing nanoslab in the clean limit, with perfectly specular 
reflection at the boundaries. Possible effects of the scat- 
tering of electrons at the interface between a substrate 
and the nanofilm will be discussed later. The transla- 
tional invariance in the perpendicular direction is broken 
due to quantum confinement and, so, we deal with a spa- 
tially nonuniform order parameter A(r) problem. As a 
consequence, the ordinary BCS self-consistent equation 
should be abandoned in favor of a more elaborate analysis 
based on the Bogoliubov-de Gennes (BdG) equations (or, 
equally, on the Gor'kov Green's function formalism). As 
shown in Ref. [16, the BdG equations are appropriate to 
describe the BCS-BEC crossover in spatially nonuniform 
fermionic systems at nearly zero temperatures. In partic- 
ular, it has been shown in Ref. (la that the BdG equations 
reproduce the Gross-Pitaevskii equation for the conden- 
sate wave function at the BEC side of the crossover. 
Therefore, if our reasoning based on _K"||/|[/| is correct, 
we can expect that the spatial profile of the wave function 
of a condensed electronic pair will change (dramatically) 
crossing over from the BCS regime (i.e., A"||/|C/| ^$> 1) to 
the BEC regime (i.e., K\\/\U\ < 1). In the BCS regime 
this pair-wave function will have many oscillations with 
the period of the Fermi wavelength A^ and will decay 
over a significantly larger distance (as compared to Xp) 
determined by the size of the extended Cooper pairs. Ap- 
proaching the BEC regime, the pair-wave function will 
have only some insignificant residual oscillations (e.g., 
due to the presence of the ultraviolet cut-off) and will 
be concentrated at short separations between electrons 
associated with the size of a molecule-like pair (it is close 
to or smaller than Xp). For a superconducting slab the 
BdG equations can be written in the form 



£i,k|ui,k) = H e \u hk ) + A|i> i]k ), 
£i,kKk) = A*|w l:k ) - H*\v it u) 



(la) 
(lb) 



where i is the quantum number associated with the 
quantum-confined motion in the z direction, k is the 
wavevector of the quasi-free motion of electrons in the 
direction parallel to the nanofilm (x and y directions); 
\ui.u) and |wi,k) are the particle-like and hole-like ket 
vectors; e^k is the quasiparticle energy. In addition, 
A = A(f), with f the electron position operator, and 



H e = H* (* stands for the complex conjugate) is the 
single-electron Hamiltonian (measured from the Fermi 
level E F ) 
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$t( r ) = 5Z [( r Kk) 7i,k,t " Kk|r) 7i,-k,; . ( 7a ) 

i,k 

$4-( r ) = X! [( r K k > 7»,-k,4. + (Ui,k|r) 7 t t kjt , (7b) 



with m e the band mass set to the free electron mass and 
V«>nf ( r ) the confining potential. Equations (fT]) are solved 
together with the self-consistency relation given by 



A ( r ) = gJ2(*\ui,k)( v iM\r), 



(3) 



where g is the coupling constant, and we limit ourselves 
to zero temperature. We ignore pairing of electrons be- 
tween different subbands, which is justified when the in- 

t2 2 

tersubband energy spacing, i.e., S ~ ^fy, * s signifi- 
cantly larger than the pairing energy (it is always true 
for ultrathin films). In this case it is possible to take 
into account only the pairing of the time reversed states, 
which means that (see, e.g., Ref. fl7l) 
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where # e |£i,k) = £i,k|fi,k), with £ i)k the single-electron 
energy measured from the Fermi level, and the factors 
Ui t k and V^k are chosen real, together with the order 
parameter. Based on Eqs. (TTJ), ^ and flU), one can find 
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with e^k = \/£ 4 2 k + Af, where the subband-dependent 

gap Aj = (fi.klAl^^k) obeys the following BCS-like self- 
consistency equation: 

A„=5> 



;,k 



' 2e 2 . 



(6) 



with $,v, t = .g/d 2 r|(r|&',k)| 2 |(r|&,k)| 2 . As usual, to 
avoid the ultraviolet divergence, the sum in Eq. ^ [and 
in Eq. ([3])] runs over the states with |£j )k | < Swd, where 
lod stands for the Debye frequency. For Pb we adopt 
huj D = 8.27 meV, see Ref.|H 

To check our arguments about the influence of the ra- 
tio Ku/U on the lateral subband-dependent size of the 
Cooper pairs, we need to study the "wave function" of 
a condensed fermionic pair, i.e., the anomalous aver- 
age of the field operators (^ (r)^»j,(r')). More precisely, 
the anomalous average is proportional to the bound-like 
eigenstate of the two-particle density matrix, which can 
be safely interpreted as the wave function of a condensed 
fermionic pair, see, e.g., Refs. [l9| andl20l This quantity 
is related to the particle-like and hole-like vectors of the 
BdG equations through the Bogoliubov canonical trans- 



where 7T and 7 are the quasiparticle (bogolon) operators. 
From Eq. £f) we find (for T = 0) 

(^t(r)^(r')) = 5> ilk V ilk (r|6,k)(^k|r'), (8) 



which can be rearranged to 
(^ t (r)4(r')> =J2^< P i(z)rt(z')r li (x-x',z-z% (9) 

i 

where <fi(z) is the single-electron wave function asso- 
ciated with the i— th QWS, and rji(x — x',y — y') — 
Vi(p) [with p = yf{x — x') 2 + (y — y') 2 } controls the de- 
cay of the fermionic-pair " wave function" in the direction 
parallel to the nanofilm. This quantity can be repre- 
sented in the form 



/c\k 
-k9{^ D -\^\) 



Jgjkp) 
K + A 2 



(10) 



with k = Jk 2 + k 2 and Jo the Bessel function of the first 
kind of order 0. Let us here introduce the QWS energy 
Xi, which makes it possible to write £i,k = ^{k 2 + k 2 ) — 
Hi, with pii = Ep — Xi. The asymptote for r]i(p) at large p 
can be calculated analytically in several interesting cases. 
When the QWS level (i.e., the subband bottom) is situ- 
ated far below the Fermi level, i.e., pi > and pi^$> A^, 
we find 

m (p)^^J (k lP )K (p/R^) [p/i?f } >l], (11) 

where Kq is the Macdonald function and R\ = Hvi/Ai, 
with Vi = hki/m the subband-dependent Fermi veloc- 
ity (for parallel motion). This is nothing more but the 
ordinary BCS behavior: first, there are fast oscillations 
with period of the subband-dependent Fermi wavelength 
Xi = 2n/ki [this comes from Jo(kip)]; second, we obtain 

the exponential overall decay governed by R\ , as seen 

from the asymptote of K (p/R\ '). 

When QWS touches the Fermi level, i.e., pi = 0, 
Eq. (|11[) does not hold any more and the behavior of r\i (p) 
changes dramatically. In this case, assuming p — > 00, we 
arrive at 



m 



Vi ( P ) ~ — Jo(p/Rn Koip/Rr') [p/Ri > 1] , 



(2)s 



(12) 

with R\ ' = h/y/mAi. As seen, the characteristic length 
controlling the decay of the pair "wave function" in 
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FIG. 1: (Color online) (a,b,c) The subband-dependent pair "wave function" r]i(p) versus p for different values of pt, i.e., 
Pi/hwD = (a), pi/fudD — 1 (b), and pi/hojo = 100 (c). Panel (d) shows how the product kiRi depends on the ratio 
Pi/hhJo (1/ki is the measure for the mean distance between electrons in a given subband and Ri stands for the pair size). The 
inset in (d) is a room around pi = 0, which shows details of a drop in the product kiRi when the bottom of the corresponding 
single-electron subband approaches the Fermi level (the BCS-BEC crossover is reached below the dashed line, i.e., kiRi < 2ir). 



Eq. (|12p is completely different from that in Eq. (ITTI) . 



? (2) 



1) 



Moreover, R\ is significantly smaller than R\ ' . As 
R\ can be well approximated by hvp/Ai, with vp — 
y / '2Ep/m, then we obtain 



rP/r 



(2) 



OC 



VEp/Ai 



(13) 



In addition, the fast oscillations present in Eq. (|lll) dis- 
appear in Eq. (|T2|) . Instead, we obtain a rather slowly 

oscillating factor Jo{p/R\ ) but the role of these oscilla- 
tions is almost negligible because they can only man- 
ifest themselves for separations of electrons in a con- 
densed pair such that r]i(p) almost approaches zero. So, 
our expectations based on the behavior of the parameter 
-ftf||/|[/| are relevant: when the subband bottom touches 
the Fermi level, the Cooper pairs associated with this 
subband nearly collapse (as compared to typical values 
of the zero-temperature BCS coherence length) in the lat- 
eral direction, shrinking by a factor of ^ Ep/ Ai ~ 10 2 . 
Such shrinking continues when p t crosses zero and be- 
comes negative. In particular, when hujp> 3> \pi\ 3> A^, 
Eq. dTO]) results in 

TH{p) ~ -^ K {p/Rf ] ) [p/Rf ] > 1] , (14) 



irh 2 



with #j 3) = 
obtain 



H/\/2m\fj,i\. As R\ ' w ^h/(2mujp))-, we 



RY ] /Rr oc y/EpJuJo/Ai 



(15) 



which can be larger than the ratio in Eq. (1131) by a fac- 
tor of y^hujj/Ai, which is about 3 — 10 for conventional 
superconductors. Thus, based on Eqs. (fl~3|) and (|15[) . we 
find a giant drop, from microns to nanometers, in the 
lateral size of a condensed fermionic pair associated with 
the single-electron subband whose bottom approaches 
Ep. This molecule- like trend in the reconstruction of 
the electronic pairing is very similar to the one predicted 
previously for high-quality superconducting nanowires. 17 



More details about r)i(p) can be obtained by numer- 
ically calculating the integral in Eq. (|10j) . Using A; = 
1.08 meV, which results from tunneling measurements of 
4 ML thick lead nanofilms^^ we investigate how the total 
profile of r\i evolves with changing pi , in addition to the 
asymptotic behavior discussed in the previous paragraph. 
Figure [T] shows numerical results for rji (p) for different 
values of pi, i.e., pi/hioo = (a), pi/hiou = 1 (b), and 
Pi/huiD = 100 (c). As seen, the numerical results fully 
support the analytical results given by Eqs. (Jill) . (fT2)l 
and (fl"4|) . In particular, we obtain the typical BCS pic- 
ture of anomalous correlations when the subband bottom 
is situated far below Ep [see panel (c)] but we approach a 
molecule-like character of the spatial distribution of elec- 
trons in a condensed pair for pi -» [see panel (a)]. Fast 
oscillations in r/i(p) at Fig. QJc) are first converted into 
slow oscillations shown in panel (b) [here Xi significantly 
increases as compared to panel (c)] and, then, are almost 
washed out in Fig.QJa). Only small residual oscillations 
can be observed in panel (a) , and the period of these os- 
cillations is about the decay radius for r\i (p) , which is in 
agreement with Eq. (fT2)l . To further highlight similari- 
ties with the ordinary BCS-BEC crossover driven by the 
Feshbach resonance in ultracold Fermi gases, we show in 
Fig. [IJd) how the product kiRi depends on pi, where 
Ri is the mean square radius calculated with the "wave 
function" r/i(p) and ki = (S^rii) 1 ' 3 , with rii the mean 
electron density in subband i. We note that here ki is 
introduced to measure the mean distance w Aj = 2ir/ki 
between the electrons in the i-th subband, notwithstand- 
ing the presence/ absence of the real Fermi motion and 
the position of the subband bottom with respect to Ep. 
When the bottom of a subband is far below the Fermi 
level this definition of ki (and of Xi) recovers the defini- 
tion used in Eq. (fTTj) . Following the results of Ref. [2l|, we 
can expect that the BCS-BEC crossover is approached 
in a given subband when 1/7T < kiRi < 2-7T. In other 
words, the size of a condensed fermionic pair Ri becomes 
smaller than the mean distance between electrons in the 
corresponding subband. Having at our disposal this in- 
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FIG. 2: (Color online) Sketch of the single-electron energy 
spectrum £;,k + Ep, with k = {k x ,k y }, for the lower (1) 
and upper (2) subbands in 4-5 ML lead single-crystalline 
nanofilnis. Two thick arrows show the subband-dependent 
Fermi wavenumbers. 



equality, we can learn from Fig.QJd), that the BCS-BEC 
crossover is definitely approached when the subband bot- 
tom approaches Ep, i.e., when \fii\/hujp) ^ 1- Here it 
is worth noting that, as seen from the inset in panel 
(d), the product kiRi slightly increases with decreas- 
ing Hi for jii < —0.2hujp). This feature appears due to 
the presence of the ultraviolet cut-off and demonstrates 
that the regime of Eq. (jT¥J) is not reached in our calcu- 
lations. The point is that Eq. (jl~4")) essentially assumes 
that Hujd ^> \im\ ^> A* but this can not be realized for 
the parameters used to calculate the data of Fig. [TJ i.e., 
for TvjJd — 8.27 meV and Aj = 1.08 meV. It is also of im- 
portance to note here that despite a significant shrinking 
in their lateral size down to Xi, the condensed pairs of 
electrons do not suffer from the Coulomb repulsive ef- 
fects. In metals such effects can be expected only for the 
separations smaller than 0.3-0.4 nm whereas in the sub- 
band whose bottom is in the vicinity of the Fermi level 
Xi ~ 5-10 nm, i.e., the mean density of electrons in the 
subband whose bottom approaches Ep is significantly 
smaller than the total mean electron density. We stress 
that this does not prevent this subband from making a 
significant contribution to the coherent phenomena, see 

secnni 

Based on the results of this section, we can conclude 
that the redistribution of the kinetic energy between 
the perpendicular and parallel spatial degrees of freedom 
leads to a significant reconstruction of the internal distri- 
bution of electrons in a condensed pair so that we find a 
clear molecule-like trend in the pairing. This behavior is 
found for the pair condensate in a subband whose bottom 
approaches the Fermi level. In total, we have a coher- 
ent mixture of different subband-dependent condensates. 
So, the question arises what are the consequences of this 
mixture in the presence of the BCS-BEC crossover-like 
behavior in one of the available subbands. This point is 
addressed in the next section. 



III. TWO CONDENSATES AND TWO 
CHARACTERISTIC LENGTHS 

In the previous section we discussed the results for the 
fcrmionic condensate associated with one of the single- 
electron subbands appearing due to quantization of the 
perpendicular electron motion. However, due to the pres- 
ence of multiple subbands we have a coherent mixture of 
subband-dependent fermionic condensates each making 
a different contribution to the total condensate and su- 
perconducting characteristics. In nanofilms, due to the 
quasi-2D character of the electron motion, such contri- 
butions are, as a rule, almost equal to each other, except 
in the case when the subband bottom becomes signifi- 
cantly higher than Ep — hwp>. In particular, when it ap- 
proaches the level Ep + hwp, , the corresponding subband 
stops to contribute due to the ultraviolet cut-off. The 
most interesting choice is realized when there are only 
two relevant subbands^ which is the case for, e.g., lead 
nanofilms with thicknesses 4 and 5 ML (see the sketch 
given in Fig. [2]). The optimal variant of tuning the sub- 
band positions is when the bottom of the upper sub- 
band is situated at Ep — hijjp,. First, the contribution 
of such a subband to the superconducting characteristics 
is almost the same as that of the lower subband with 
the bottom far below Ep. Second, the lateral size of 
the condensed fcrmionic pairs associated with the upper 
subband in this case only slightly increases as compared 
to Hi — 0, see Fig. E^d). Our numerical results for this 
variant are shown in Fig. [3] for lead nanofilms with thick- 
nesses 4 ML (the left panel) and 5 ML (the right panel) . 
1 ML corresponds to 0.286 nm in our calculations^ Here 
the quantity 



w{R) = 1 



ffd 3 rdV(^(r)^(r')) 

P<R 

jJd3rdV(^ t (r)^(r')> ' 



(16) 



i.e., the probability of finding a Cooper pair with the lat- 
eral size larger than R, is plotted versus R. In both 
panels of Fig. [3] the different curves are for: (1), the 
results for the upper subband, whose bottom touches 
Ep — hup, [when keeping only its contribution to Eq. (J9J] ; 
(2), the total results of both contributing subbands; and 
(3), the data for the lower single-electron subband whose 
bottom is far below Ep. Here the coupling constant g is 
chosen such that the experimental critical temperature 
T c , as reported in Ref.[3, is obtained. For illustrative pur- 
poses, the confining interaction was taken as zero inside 
the nanofilm and infinite outside. The Fermi level Ep is 
adjusted so that to get the second QWS at Ep — ftuin, 
i.e., Ep — 3<2 + fkup- As seen from Fig. 02 the results 
for 4 ML and 5 ML are quite close to each other. If we 
adopt that the lateral diameter of the fermionic pair is 
defined by w(R) = 0.32 (recall that for the normal distri- 
bution 1— f_° dTp„ orm (T) = 0.32), then we find that the 
Cooper-pair radius associated with the lower subband is 
Riow = 72 nm for 4 ML and 62 nm for 5 ML. The minor 
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FIG. 3: (Color online) The probability of finding a Cooper 
pair with lateral size larger than 7? for Pb superconducting 
nanofilms with thickness 4 ML (the left panel) and 5 ML (the 
right panel). The results for a condensed fermionic pair in 
the upper subband are given by curve 1; the total probability 
for the both contributing subbands is represented by curve 
2; the results for the Cooper pairs associated with the lower 
subband are given by 3. The Fermi level is adjusted so that 
the bottom of the upper of the two available subbands is at the 
energy Ef — hcun ■ The dotted line in both panels correspond 
to w = 0.32. 



difference here is due to slightly different Fermi veloci- 
ties (in the parallel direction) and critical temperatures. 
As for the pair radius for the upper subband, we obtain 
Rup = 7.6 nm for (a) and R up = 7.8 nm for (b). These 
results are well explained by the trend given by Eq. (|12[) , 
i.e., R up ex h/y/mk B T c , with T c = 6.7 K for 4 ML and 
6.3 K for 5 ML (see Ref. y|). For the total radius we have 
for curve (2) R to t = 34 nm for panel (a) and i? t ot = 32 nm 
for panel (b) . The drop of i? t ot as compared to the bulk- 
like result for the lower subband is due to the contri- 
bution of the extremely small pairs associated with the 
upper subband. It is worth noting that the introduction 
of the characteristic pair size for a coherent mixture of 
two different condensates is somewhat conventional. As 
seen from Fig. [3J curve (2) exhibits a bimodal behavior 
with a clear crossover from the short-range regime gov- 
erned by the upper band to the long-range one controlled 
by the lower band. Such a behavior can not be, of course, 
accurately specified by one characteristic length. This is 
in agreement with the recent analytical results for the 
extended Ginzburg-Landau formalism for multiband su- 
perconductors 11 and with numerical investigations of the 
coherent properties of MgB2, see, e.g., Ref. |2J. Thus, 
R to t measures the short-to- long crossover in w(R) rather 
than the size of a condensed fermionic pair. Neverthe- 
less, the probability to pick up a condensed pair with 
lateral size smaller than R tot is about 0.68. Moreover, as 
w(R) decreases relatively slow around -Rtot, the probabil- 
ity to find a condensed pair with size smaller than 10 nm 
is about 0.5: almost half of the condensed fermionic pairs 
are smaller than 10 nm in lateral size. 

What are possible experimental evidences of such 
shrinking of the Cooper pairs due to the BCS-BEC 
crossover-like behavior in single-crystalline nanofilms? 
First of all, this is a competition of different characteristic 



spatial lengths of two condensates. It was recently shown 
that such a competition can lead to the appearance of 
long-range attraction between vortices in a two-band su- 
perconductor, see Ref. [25). In this case the Abrikosov lat- 
tice of vortices formed in the perpendicular magnetic field 
will melt, with possible formation of stripes and clusters 
of vorticesJ^ Second, the profile of the order parameter 
in the core of the vortex will be sensitive to the presence 
of the condensate with extremely small spatial coherence 
length, which can be probed via the scanning tunneling 
microscopy, see, e.g., Ref. [2- Third, we can also expect 
significant proliferation of pair fluctuations in the pres- 
ence of a single-electron subband with the bottom next 
to Ep. Such fluctuations will result in the suppression 
of the DOS at the Fermi level above the critical tem- 
perature^, i.e., the pseudogap behavior. However, it 
is difficult to say whether or not these pairing fluctua- 
tions make a contribution to the pseudogap-like behav- 
ior revealed recently in single-crystalline lead nanofilms 6 : 
other mechanisms can also be involved (see discussion in 
Ref. Mj and, so, a more involved analysis is needed. In 
addition, Josephson physics and Andreev states will be 
also very sensitive to the appearance of the molecule- 
like trend in pairing and to the competition of different 
length-scales. In particular, these important things can 
be probed with hybrid superconducting devices similar 
to, e. g., a carbon nanotube joining a superconducting 
loopi^s In our case the single-crystalline metallic nanofilm 
is supposed to transmit the current between supercon- 
ductors over mesoscopic distances (to keep the nanofilm 
in the normal state it is necessary to have the critical 
temperature of leads higher than that of the nanofilm) . 



IV. CONCLUSION 

In conclusion, we demonstrated that quantum confine- 
ment of the perpendicular motion of electrons in single- 
crystalline metallic superconducting nanofilms leads to 
the formation of a nanoscale multi-band(subband) super- 
conductor. The position in energy of each subband can 
vary significantly depending on the fabrication process, 
and it is possible that the bottom of one of such sub- 
bands is situated in the vicinity of the Fermi level. We 
showed that the character of the superconducting pair- 
ing in such a subband changes dramatically and exhibits 
a clear molecule-like trend, which is very similar to the 
BCS-BEC crossover but now driven by the perpendicular 
size-quantization. 

Though we have considered the system in the clean 
limit, our results will also be relevant in the presence 
of a moderate disorder. The main issue for disorder in 
high-quality metallic superconducting nanofilms comes 
from the scattering of electrons at the interface between 
a substrate and the nanofilm, which leads to a broad- 
ening of the QWS levels. Though this will smoothen a 
drop in the lateral size of the Cooper pairs, our results 
given in Fig. [^d) make it possible to expect that the ef- 



feet in question will survive to a great extent until the 
broadening essentially exceeds bfrwu ~ 40 meV. Here it 
is worth noting that in experiments with lead nanofilms, 
the level broadening of QWS can be reduced, in the pres- 
ence of a crystalline interface, down to 30-40 meV (at zero 
temperature and at energies close to Ef), see Ref. 1271 
This does not exceed 5 hujo mentioned above, which gives 
quite an optimistic view on the prospects of an experi- 
mental observation of the effect in question. We would 
also like to note that the line width for the perpendic- 
ular discrete levels increases with temperature due to 
electron-phonon scattering. Though such an increase is 
not very pronounced for temperatures around T c , i.e., 
it is of about 20meV, see Ref. [271 this will cause an 
additional temperature-dependent smoothing. So, the 
optimal choice to probe the BCS-BEC-like crossover in 
single-crystalline nanofilms is to stay at nearly zero tem- 
peratures. 

We would also like to highlight strong similarities be- 
tween our findings and recent angle resolved photoe- 
mission results 13 for multiband superconducting iron- 
chalcogenide FeSe x Tci_ x . It was found in Ref. |T3| that 
the investigated sample of FeSe x Tei_ x does not lay deep 
in either the BCS or BEC regimes but lays instead in the 



BCS-BEC crossover domain. As demonstrated above, 
the same conclusion holds for 4 ML and 5 ML thick lead 
nanofilms when the bottom of one of the relevant single- 
electron subbands approaches Ep. Thus, one can in gen- 
eral conclude that the multiband superconductors (both, 
bulk and with quantum-confined induced subbands, see 
also Refs . KL7I and \2&i) offer new important possibilities to- 
wards further understanding the physics of the BCS-BEC 
crossover in solid-state materials. 
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